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ABSTRACT 

Context. This paper is the third in a series of investigations by the authors. 

Aims. The nature of fast magnetoacoustic and Alfven waves is investigated in a 2D = plasma in the neighbourhood of two dipoies. 
Methods. We use both numerical simulations (two-step Lax-Wendroff scheme) and analytical techniques (WKB approximation). 
Results. It is found that the propagation of the linear fast wave is dictated by the Alfven speed profile and that close to the null, the 
wave is attracted to the neutral point. However, it is also found that in this magnetic configuration some of the wave can escape the 
refraction eff'ect; this had not been seen in previous investigations by the authors. The wave split occurs near the regions of very high 
Alfven speed (found near the loci of the two dipoies). Also, for the set-up investigated it was found that 40% of the wave energy 
accumulates at the null. Ohmic dissipation will then extract the wave energy at this point. The Alfven wave behaves in a diff'erent 
manner in that part of the wave accumulates along the separatrices and part escapes. Hence, the current density will accumulate at 
this part of the topology and this is where wave heating will occur. 

Conclusions. The phenomenon of wave accumulation at a specific place is a feature of both wave types, as is the result that a fraction 
of the wave can now escape the numerical box when propagating in this magnetic configuration. 
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1. Introduction 

The solar corona is a complicated, inhomogenous magnetised 
environment. It is also know that there are MHD wave mo- 
tions in the corona and that the environment plays a key role in 
their propagation. Wave motions have been observed with the 
SOHO and TRACE satellites, for example slow MHD waves 
(Berghmans & Clette, 1999; De Moortel et al, 2000) and fast 
MHD waves (Nakariakov et ah, 1999). To begin to understand 
this magnetic environment, it is useful to look at the structure 
■ (topology) of the magnetic field itself. Potential field extrapo- 
" lations of the coronal magnetic field can be made from photo- 
spheric magnetograms. Such extrapolations show the existence 
' of an important feature of the topology; null points. Null points 
' are points in the field where the Alfven speed is zero. Detailed 
investigations of the coronal magne tic field, using such pote ntial 
field calculations, can be found in |Beveridge et al. (2002)) and 
[Brown & Priest (200 1) [ 

McLaughlin & Hood (2004) (hereafter referred to as Paper I) 
found that for a single 2D null point, the fast magnetoacoustic 
wave was attracted to the null and the wave energy accumulated 
there. In addition, they found that the Alfven wave energy ac- 
cumulated along the separatrices of the topology. McLaughlin 
& Hood (2005) (hereafter Paper II) found that the key results 
from Paper I carried over from the simple 2D single magnetic 
null point to the more complicated magnetic configuration of 
two null points. The aim of this paper (eflfectively the third in a 
series) is to see if their ideas carry through to a more realistic 
magnetic configuration (that still involves a magnetic null). 
Waves in the neighbourhood of a single 2D null point have been 
investigated by various authors. [Bulanov & Syrovatskii (1980)| 
provided a detailed discussion of the propagation of fast and 
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Alfven waves using cylindrical symmetry. In their paper, har- 
monic fast waves are generated and these propagate towards the 
null point. However, the assumed cylindrical symmetry means 
that the disturbances can only propagate either towards or away 
from the null point. |Craig & Watson (1992)| mainly consider the 
radial propagation of the m = mode (where m is the azimuthal 
wavenumber) using a mixture of analytical and numerical so- 
lutions. In their investigation, the outer radial boundary is held 
fixed so that any outgoing waves will be reflected back towards 
the null point. This means that all the energy in the wave motions 
is contained within a fixed region. They show that the propaga- 
tion of the m = wave towards the null point generates an expo- 
nentially large increase in the current density and that magnetic 
resistivity dissipates this current in a time related to log rj. Their 
initial disturbance is given as a function of radius. In this paper, 
we are interested in generating the disturbances at the boundary 
rather than internally. Craig and McClymont (1991, 1993) and 
Hassam (1992) investigate the normal mode solutions for both 
m = and m ^ modes with resistivity included. Again they 
emphasise that the current builds up as the inverse square of the 
radial distance from the null point. However, attention was re- 
stricted to a circular reflecting boundary. 

The propagation of fast magnetoacoustic waves in an in- 
homogeneous coronal plasma has been investigated by 
Nak ariakov & Roberts (1995)[ who showed how the waves are 
refracted into regions of low Alfven speed. In the case of null 
points, it is the aim of this paper to see how this refraction pro- 
ceeds when the Alfven speed actually drops to zero. 

The propagation of fast waves in 2D coronal arcades has been in- 
vestigated by [Oliver g/. (1993)| When the background Alfven 
speed is either constant or decreases monotonically with height, 
the waves are leaky ( D iaz et al. 2 005 ). However, when the den- 
sity is chosen so that there exists a local minimum in the Alfven 
speed, the fast waves can be trapped ([Verwichte et al. 20051 
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|Brady & Arber 2Q()5| ). The present paper investigates the nature 
of fast wave propagation when there is a null point at which the 
Alfven speed actually drops to zero at a local minimum. The 
same idea of trapping at least a fraction of the wave is found. 
Thus, fast waves in 2D coronal arcades can be trapped by either 
a density enhancement or the presence of a null point. 
The paper has the following outline. In Section[2l the basic equa- 
tions are described. The results for an uncoupled fast magne- 
toacoustic wave are presented in Section [3l This section dis- 
cusses fast wave propagation with a pulse generated on the lower 
boundary, and numerical and anaytical results for three scenar- 
ios are presented. Section [4] discusses the propagation of Alfven 
waves and the conclusions are given in Section!?] 



2. Basic Equations 

The usual MHD equations for an ideal, zero p plasma appropri- 
ate to the solar corona are used. Hence, 



idy \ 1 

p — -h(v-V)v = -(VxB)xB, 

\dt I /I 

5 = Vx(vxB) , 

ot 

^ + V . (pv) = , 
ot 



(1) 
(2) 
(3) 



where p is the mass density, v is the plasma velocity, B the 
magnetic induction (usually called the magnetic field) and yu = 
471 X 10"'^Hm"^ the magnetic permeability. The gas pressure and 
the adiabatic energy equation are neglected in the zero JS approx- 
imation. The magnetic diff'usivity is neglected in the ideal ap- 
proximation. 



2.1. Basic equilibrium 

Paper II showed that the key results from Paper I carried over 
from the simple 2D single magnetic null point to the more com- 
plicated magnetic configuration of two null points. However, the 
choice of magnetic fields used in both these papers are only valid 
in the neighbourhood of the null point(s), because as x and z get 
very large, Bq also gets large, which is unphysical. Thus, in this 
paper we now want to investigate the behaviour of MHD waves 
near a more realistic magnetic configuration (and one that still 
involves a magnetic null). One such choice is a magnetic field 
created by two dipoles, of the form Bq = B{)L^ (Bx, 0, 3^), where: 



-{x + a) . -{x- a) + 

[(x -h a) 

2(x-\-a)z 



] [(■' 



■af^z^ 



2(x- a)z 



[(jc + a)2 + z2] [(jc-a)2+z2] 



(4) 



where is a characteristic field strength, L is the length scale 
for magnetic field variations and 2a is the separation of the 
dipoles. We choose a = 0.5 L in our simulations. This magnetic 
field can be seen in Figure [T] It comprises of four separatrices 
and an X-point located at (x, z) = (0, a). Note that as x or z gets 
very large, the field strength becomes small (this is a more phys- 
ical field than those previously investigated in Papers I and II). 



2.2. Linearised equations 

As in Papers I and II, the linearised MHD equations are used to 
study the nature of wave propagation near the two dipoles. Using 
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Fig. 1. Our choice of equilibrium magnetic field. 



subscripts of for equilibrium quantities and 1 for perturbed 
quantities, the linearised equation of motion becomes 



Po- 



VxB 



xBn 



the linearised induction equation 
^Bi 

^ = Vx(vixBo), 
and the linearised equation of mass continuity 

^ + V.(poVi) = 0. 
ot 



(5) 



(6) 



(7) 



We will not discuss equation d?]) further as it can be solved 
once we know vi. In fact, it has no influence on the mo- 
mentum equation (in the zero j3 approximation) and so in ef- 
fect the plasma is arbitrarily compressible (Craig & Watson 
1992). We assume the background gas density is uniform 
and label it as pp. A spatial variation in po can cause 
phase mixing ([Heyvaerts & Priest 1983llDe Moortel et al. 19991 
IHood etal 20021) : 

We now consider a change of scale to non-dimensionalise; let 
vi = vv;. Bo = 5oB;, Bi = 5oB;, X = Lx\ z = Lz\ V = iV* 
and t = tt*, where we let * denote a dimensionless quantity and 
Y, Bq, L and t are constants with the dimensions of the variable 
they are scaling. We then set = v and y = L/t (this sets v as 



a sort of constant background Alfven speed). This process non- 
dimensionalises equations ([5]) and Q, and under these scalings, 
t* = I (for example) refers to ^ = ? = L/v; i.e. the time taken to 
travel a distance L at the background Alfven speed. For the rest 
of this paper, we drop the star indices; the fact that they are now 
non-dimensionalised is understood. 

The ideal linearised MHD equations naturally decouple into two 
equations for the fast MHD wave and the Alfven wave. The slow 
MHD wave is absent in this limit and there is no velocity com- 
ponent along the background magnetic field, as can be seen by 
taking the scalar product of equation ([5]) with Bq. 
The linearised equations for the fast magnetoacoustic wave are: 



db, dbx 
— = v,(x,z)\- — 



dV 

~dt 

dt 



dx 



dV 



dt 



dz 
dV 
dx 



(8) 
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where Bi = (bx,0,b^), the Alfven speed va(x,z) is equal to 
^Bl + and the variable V is related to the perpendicular ve- 
locity; V = [(vi X Bo) • Cy]. These equations can be combined to 
form a single wave equation: 



(9) 



The linearised equations for the Alfven wave, with vi = (o, Vy, 6j 
and Bi = (o, by, o) are: 



5v„ 



db. 



+ B, 



dt dx ' ^ dz ' dt 

which can be combined to form a single wave equation: 



dVy dVy 

ox oz 



dfi 



ox oz 



(10) 



(11) 



It is worth noting that the zero JS approximation will be in- 
valid near the null points, where the pressure terms become 
important. In fact, a finite JS model would introduce slow 
waves, which would be d riven by the excited fast waves 
( [McLaughlin & Hood 2006| ). Also, the model is linear but the 
inhomogenity of the medium would le ad to nonlinear coupling 
of the modes ( [Nakariakov a/. , 1997| ). 

2.3. Alfven speed profile 

We shall be solving equation (|9]) numerically in the follow- 
ing sections and sending a planar wave into the system from 
the lower boundary of our numerical box. However, due to our 
choice of equilibrium magnetic field, there are two singularities 
in the field at (x, z) = (±a, 0) and hence ^ oo at these points. 
If we were to generate a wave along z = 0, our numerical scheme 
would be unable to deal with this extreme speed diff'erential. 
Thus, in this chapter we will generate our waves not along z = 
but along z = Zo, where zo is small (and so close to the z = line). 
This choice still starts the waves in a region of varying Alfven 
speed, but reduces the massive diff'erential to something easier to 
work with numerically. Thus, starting the waves at z = Zo results 
in very little loss of information about the system. In addition, 
the CFL condition dictates that our timestep has to be very small 
(because our speed is so large in some places). Thus, to reduce 
the computing time and with very little loss of understanding, 
Zo = 0.1 was set in the experiments (recall a = 0.5 in our simu- 
lations). 

Figure [21 shows a shaded surface (left) and colour contour (right) 
of v\(x, z) = bI(x, z)-^B^(x, z). The shaded surface clearly shows 
that the velocity profile changes massively across the magnetic 
region and reaches a maximum at (x, z) = (±a,0). The contour 
(right) is colour coded; red represents values < < 1, blue 
represents 1 < < 20 and black represents v| > 20. Thus, 
we can see that around the X-point there is a small island of 
low Alfven speed. The white curve shows the skeleton of the 
magnetic field (equation [4]). 



3. Fast waves 

We now look at the behaviour of fast MHD waves in the neigh- 
bourhood of two dipoles. There are three scenarios which we 
shall visit in order to fully understand the system. Firstly, the 
two dipoles create a natural X-point and we want to know if the 



fast wave wraps around this null point in a similar way to that de- 
scribed in Paper I. Thus, Simulation One describes a localised, 
linear fast wave approaching this X-point. This magnetic con- 
figuration also contains regions of very high Alfven speed along 
the lower boundary. Thus, Simulation Two describes a linear fast 
wave generated in this area, but localised to one side of the (sym- 
metric) system. Simulation Three describes a scenario similar to 
Simulation Two, but now with a pulse generated along the entire 
lower boundary. Simulations One and Two will build towards an 
understanding of the complete system of Simulation Three. 
All three simulations generate a linear fast wave pulse on the 
lower boundary, which then propagates up into the magnetic 
configuration. The lower boundary was set at zo in all three ex- 
periments, where zo = 0.1, a = 0.5, oj = and we focus our 
attention in the region -1 < x < 1 and zo < z <2. 

3.1. Simulation One 

We solve the linearised MHD equations for the fast wave (equa- 
tions [8]), numerically using a two-step Lax-Wendroff' scheme. 
The numerical scheme is run in a box with -0.6 < x < 0.6 
and Zo < z < I, using the magnetic field shown in Figure [TJ We 
initially consider a single wave pulse coming in from the bottom 
boundary, localised along -0.2 < x < 0.2, z = Zo = OA. The 
boundary conditions were set such that: 



V(x,zoJ) = 



smojt sin (x -h l)j for 



x^-0.6 





dx 



0.2 < X < 0.2 

0<t<^ 
otherwise 



x^O.6 



0. 



Tests show that the central behaviour is largely unaff'ected by 
these choices of side and top boundary conditions. The other 
boundary conditions on the perturbed magnetic field follow from 
the remaining equations and the solenodial condition. 
We find that the linear, fast magnetoacoustic wave propagates up 
from the lower boundary (from zo) and begins to wrap around the 
null point. This can be seen in Figure [3l In the early subfigures, 
we can see that the pulse travels faster at the ends (wings) than 
in its centre, due to the variation in Alfven speed. The wave then 
continues to refract around the X-point, and wraps itself around 
it. The Alfven speed is zero at (x, z) = (0, 0.5) (the null) and so 
the fast wave cannot cross this point. Consequently, the null acts 
as an immovable barrier to the propagation of the wave and thus 
is a focus for the wrapping around eff'ect. The wave continues 
to wrap around the null and becomes thinner. The wave (even- 
tually) accumulates very close to the X-point. This is a similar 
eff'ect to that described in Paper I. This refraction eff'ect is a key 
feature of fast wave propagation. 

Since the Alfven speed drops to zero at the null point, the wave 
never reaches there, but the length scales (this can be thought 
of as the distance between the leading and trailing edges of the 
wave pulse) rapidly decrease, indicating that the current (and all 
other gradients) will increase. 

3.2. Simulation Two 

In this second simulation, we again solve the linearised MHD 
equations for the fast wave numerically (using the same two-step 
Lax-Wendroff' scheme). We again consider a single wave pulse 
coming in from the lower boundary, but now localised along 
-l<x<0, z = Zo = 0.1. This allows us to explore a diff'er- 
ent topological area than that of Simulation One. The numerical 
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Fig. 2. Left shows a shaded surface of v^. Right shows a colour contour of v^; red represents values < < 1, blue 1 < < 20 
and black v\ > 20. The white curve shows the magnetic skeleton. Both representations are cut-off at z = zo = 0.1. 
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Fig. 3. Contours of numerical simulation V for a fast wave sent in from lower boundary for -0.2 < x < 0.2, z = OA and its resultant 
propagation at times (a) ^=0.05, (b) t=OA, (c) t=0A5, (d) t =0.2, (e) t=025 and (/) ^=0.3, (g) t=035, (h) t=OA, (/) ^=0.45, (7) 
t =0.5, (k) ^=0.55 and (/) t=0.6, labelling from top left to bottom right. The black curve shows the magnetic skeleton with the null 
point located at (x, z) = (0, 0.5). 



scheme is run in a box with -1.4 < x < 1.4 and zo < z < l.l. 
The boundary conditions were set such that: 



V(x,zoJ) = 




The results can be seen in Figure IH We find that the localised 
pulse starting between -1 < x < rises, but not all parts rise at 
the same speed. The central part of the pulse rises much faster, 
with the maximum occuring over x = -0.5. This is due to the 
high Alfven speed localised in that area (as seen in Section [23]) . 
This deforms the wave pulse from its original planar form. This 
area of high speed is a new eff'ect not previously seen in Papers 
I or II. 



Part of the wave pulse also approaches the null (extreme right of 
the generated wave), and this part appears to get caught around 
the X-point. As in Simulation One, the wave cannot cross the 
null and refracts around it. This also occurs here; part of the 
wave gets caught and wraps around the null (again and again). 
However, since the wave is being stretched out and thinned as 
it encircles the null, our resolution runs out at after some point. 
Thus, we cannot clearly see the wave wrapping around the null in 
Figures m However, in Figure [5] we can see a blow-up of the re- 
gion around the X-point and thus can see that the eff'ect proceeds 
in a similar way to the asymmetric pulse described in Paper I. 
We have plotted (x, z) (-x, -z) to help the comparison (i.e. 
reflected the image in the line z = -x). Thus, this part of the 
fast wave pulse accumulates at the X-point. The resolution in 
Figure [5]is coarse since it is a blow-up from the resolution of the 
simulation. 
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Fig. 5. Contour of the numerical simulation of V for a fast wave 
sent in from lower boundary for -1 < x < after time t = 0.44. 
Figure shows blow-up of region around the null point. The image 
has been reflected in the line z = -x. 



Meanwhile, the rest of the wave continues to propagate upwards 
and spread out. This part proceeds at a slower speed than be- 
fore, again due to the change in Alfven speed. Once above the 
magnetic skeleton, the wave appears to simply continue to rise 
and spread out. The wave is no longer heavily influenced by the 
magnetic null point; it has, in efl'ect, escaped the refraction ef- 
fect. This is again a new feature not seen previously in Papers I 
or II. 

However, since part of the wave is wrapping around the null and 
a second part is rising away from the magnetic region, there must 
come a point where the wave splits. This can be seen in the lower 
set of subfigures in Figure IH Turning our attention to the be- 
haviour on the right of the null, we see that the wave continues 
to spread out; part of it appears to be travelling to the top right 
corner and part of it is hooked around the null. We see the wave 
is stretched between these two goals. The wave then splits, with 
part of it ultimately wrapping around the null (and continuing 
much in the same way as the asymmetric wave pulse in Paper I) 
and the other part of the wave then continues to spread out and 
travel away from the magnetic region. This can be seen in the 
bottom two rows of subfigures. The actual wave splitting itself 
occurs when the wave enters the right hand region of high Alfven 
speed (large v\ around (±0.5, 0)). This thins the wave massively 
and forces the two parts to split. 

Thus, a wave pulse starting between -1 < x < travels into the 
magnetic region (unevenly due to the change in speed profile). 
Part of the wave experiences a slight refraction eflfect due to the 
Alfven speed profile, but ultimately spreads out and propagates 
away. However, part of the wave is caught by the null and re- 
fracted into the X-point. This part of the wave accumulates at 
the X-point. 



3.2.1 . WKB approximation 

It is clear from Simulation Two that quite a complicated split- 
ting eff'ect is in operation. We can, however, use a WKB solution 
to gain extra insight into the numerical simulation. Substituting 
y _ ^i(p(x,z) . ^-lojt jj^^Q equation (|9]) and assuming 6l) » 1 
(WKB approximation) leads to a first order PDF of the form 
*F (x, z, (p, p, q) - 0. Applying the method of characteristics, we 
generate the equations: 



J0 

ds 
dp 
ds 



I d 



dx 
ds 
dz 
ds 



(12) 



where p = q = Bx and are the components of our 
equilibrium field, oj is the frequency of our wave and s is some 
parameter along the characteristic. These five ODFs were solved 
numerically using a fourth-order Runge-Kutta method. 
Thus, by only plotting the fieldlines coming from -1 < < 0, 
we can construct the analytical equivalent of Simulation Two. 
This can be seen in Figure [6l Also, with correct choices of s 
the solution can be directly compared to our numerical solution. 
The agreement between the numerical work and WKB approxi- 
mation is very good (the contours essentially lie on top of each 
other). Following the splitting of the fast wave in Figure [6]is par- 
ticularly interesting. In the figure, the position of each wavefront 
is indicated by many tiny crosses. This allows the reader to more 
clearly follow see the wave stretching and splitting. 
We can also use our WKB solution to plot the particle paths of 
individual elements from the initial wave. This can be seen in 
Figure [71 

In the left subfigure, we see the particle paths for starting points 
of -1 < xq < set at intervals of 0.01. The line xq = -0.4 
has been removed (this starting point leads to a spiralling par- 
ticle path that would confuse the plot). We see that the lines 
for xq < -0.5 do not appear to be influenced by the null and 
simply propagate at varying angles. The particle paths for start- 
ing points of xq > -0.5 and (approximately) xq < -0.4 are 
influenced by the null, but only in so much as to deflect the 
ray. However, for starting points greater than (approximately) 
Xq = -0.4, the particle paths spiral in towards the X-point and 
are trapped. Thus, there is a critical starting point that divides 
these two types of behaviour. The central figure shows the par- 
ticle paths for xq = -0.403 (red) and xq = -0.402 (blue). We 
can see that the critical starting point is thus somewhere in the 
middle of these two. Increasing the resolution shows the criti- 
cal starting point to be xq = -0.40236 to five decimal places. 
Finally, the right hand figure shows the particle path for a start- 
ing point of Xq = -0.4 (note the change in axes). This clearly 
demonstrates the particle spiralling into the null. 

3.3. Simulation Three 

Our final simulation again solves equations ([5]) using our two- 
step Lax-Wendroff' scheme, but this time for a wave pulse com- 
ing in from the entire lower boundary. Using the magnetic field 
in Figure [TJ the boundary conditions were set such that: 



V(x,zoJ) 

dv 

dx 



-l<x<l 

[2 ' ] < t < ^ 

~ ~ OJ 

otherwise 



sin ojt sin [f (-^ + 1)] for 



= 0, 



x=-1.4 



dv 

dx 



= 0, 



x=lA 



dv 



= 



z=1.7 



We find that the linear, fast MHD wave travels up from the lower 
boundary and has its propagation influenced by the magnetic 
configuration. The results can be seen in Figure [H The resultant 
propagation is identical to the behaviour described in Simulation 
Two, since the magnetic region is symmetrical about x = 0. 
Initially, we see the wave distorts due to the two regions of en- 
hanced Alfven speed. Then, the central part of the wave slows 
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Fig. 4. Contours of numerical simulation of V for a fast wave sent in from lower boundary for -1 <x<0, z = 0.1 and its resultant 
propagation at times {a) ^=0.055, (b) ^=0.11, (c) ^=0.165, {d) t =0.22, {e) t=0.215 and (/) t=033, (g) r=0.385, (h) t=0A4, (/) 
^=0.495, 0') t =0.55, (k) ^=0.605, (/) ^=0.66, (m) ^=0.72, (n) ^=0.74, (o) ^=0.76, (p) t =0.78, (q) r=0.8 and (r) ^=0.82, (s) t=OM, (t) 
t=OM, (u) r=0.88, (v) t =0.9, (w) ^=0.92 and (x) ^=0.94, labelling from top left to bottom right. 



down and refracts around the null point (as seen in Simulation 
One). The rest of the wave (the wings) continue to propagate 
upwards and because the wave is spreading out (as seen most 
clearly in the particle paths of Figure [7]), the wave gives the im- 
pression of pivoting about the null. As this pivoting continues, 
the two wings eventually cross each other. However, due to the 
linear nature of the model, both wings pass each other without 
interacting. 

The two wings continue to pivot, and near to the null part of the 
wave wraps tightly around it. This can be seen clearly in Figure 
[9j This shows a blow-up of the wave accumulation around the 
null. Even though the resolution is coarse in this area, it is clear 



that this is the same behaviour seen in Paper I for the fast wave 
wrapping around the simple X-point. In Figure [9] we have let 
(x, z) (-X, -z) to help the comparison (i.e. again reflected the 
image in the line z = -x). 



In the lower half of subfigures of Figure[8l the wave is eventually 
stretched so much in the region -0.6 < x < 0.6, 0.2 < z < 0.5 
that it splits; part of the wave remains close to the X-point (and 
ultimately accumulates there) and part of the wave is repelled 
by the regions of high Alfven speed and is ejected away from 
the magnetic skeleton. These parts then continue to propagate 
upwards and outwards, forming a distinctive cross shape. 
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Fig. 6. Analytical solution of V for WKB approximation of a fast wave sent in from lower boundary for -1 < x < 0, z = 0.1 and its 
resultant propagation at times (a) ^=0.055, (b) t=0.n, (c) t=0A65, (d) t =0.22, (e) t=0.215 and (/) ^=0.33, (g) r=0.385, (h) t=0A4, 
(0 ^=0.495, (7) t =0.55, (k) ^=0.605, (/) ^=0.66, (m) t=0J2, (n) ^=0.74, (o) r=0.76, (p) t =0.78, (q) r=0.8, (r) r=0.82, (s) ^=0.84, (0 
r=0.86, (u) r=0.88, (v) t =0.9, (w) t=0.92 and (x) ^=0.94, labelling from top left to bottom right. The curves consisting of the crosses 
represent the front, middle and back edges of the WKB wave solution, where the pulse enters from the top of the box. The magnetic 
skeleton is also shown. 



3.3.1 . WKB approximation 

We can again use the WKB approximation to find a semi- 
analytical solution, and thus compare it to Simulation Three. 
In Section 13.2.11 we found a WKB approximation for our fast 
wave equation (equations [12]). Thus, by plotting the fieldlines 
coming from -1 < xq < 1, we can make a direct comparison 
with Simulation Three. This can be seen in Figure [TOl where the 



wavefronts are plotted at the same times as Figure [51 Again, the 
agreement is very good. As before, each wavefront consists of 
many tiny crosses. 

Note that in the numerical simulation, the fast wave diff'uses 
slightly and resolution runs out, hence some parts of the wave 
are unresolved fully. This is not the case in the WKB approxi- 
mation, where the wavefronts are crisp. Thus, the numerical and 





-0.10 -0.05 0.00 0.05 0.10 

X 

Fig. 9. Contour of the numerical simulation of V for a fast wave 
sent in from lower boundary for -1 < x < 1 after time t = 0.44. 
Figure shows blow-up of region around the null point. The image 
has been reflected in the line z = -x. 



WKB results show the same behaviour, but have small difl'er- 
ences close to the leading and trailing wavefronts. However, the 
comparison with the middle WKB wavefront is excellent. 
We can also use our WKB solution to map out the particle paths 
of the wavefront. This can be seen in Figure[TTl Here, we can see 
particle paths for starting points of - 1 < xq < 1 set at intervals of 
0.01. The blue lines shows starting points -0.4 < x < 0.4, green 
lines show starting points of xq = ±0.4, and all other starting 
points (xq < -0.4, xo > 0.4) are in red. We can clearly see that 
the critical points (where the wave splits to head to the null or 
away from the magnetic skeleton) are around xq = ±0.4. 
Hence, we see that for fast magnetoacoustic waves in the neigh- 
bourhood of two dipoles, a fraction of the wave (and thus wave 
energy) is trapped by the null with the rest escaping. The fraction 
captured by the null will be of the form j^x^^if (a, zo), where L is 
the length of our lower boundary (-L < x < L) and x^^/^ is the 
critical starting point that divides the particle paths into those 
that spiral into the null and those that do not (e.g. the starting 
point connected to the green lines in Figure [TT]). For the linear 
fast wave investigated in this paper; L = I, z = Zo = 0.1, a = 0.5 
and x^^/^ = 0.4, which means 40% of the wave is trapped by 
the null. Hence, 40% of the wave energy accumulates at the null 
and ohmic dissipation will extract the energy in the wave at this 
point. 

4. Alfven waves 

We now turn our attention to studying the behaviour of the 
Alfven wave in the neighbourhood of two dipoles. Again, we 



2.0 




X 



Fig. 11. Particle paths for starting points of -1 < < 1 set at 
intervals of 0.01. Blue lines shows starting points -0.4 < x < 
0.4. Green lines show starting points of xq = ±0.4. All other 
starting points (xq < 0.4, xq > 0.4) are in red. 



will be launching a wave into the system from the lower bound- 
ary (z = zo). However, we already have a good idea of how the 
wave will behave due to the work carried out in Paper II, i.e. that 
the linear Alfven wave is always confined to propagate along the 
field line it starts on. 

Hence, we will not carry out three similar simulations (as we 
did above for the fast wave). Instead, we will look a one case 
for the Alfven wave and find that the behaviour is indead of the 
same nature as found in Paper II. Thus, this section is included 
for completeness, but no further cases are shown here (since any 
other experiments would just show the same type of predictable 
behaviour). 

4.1. Numerical simulation 

The equations describing the behaviour of the Alfven wave 
(equations [TOl) were solved numerically using a two-step Lax- 
Wendroff' scheme. We initially consider a box (-1 < x < 1, 
Zo < z <2) using the magnetic field shown in Figure [T] We con- 
sider a single wave pulse coming in from the bottom boundary, 
localised along -1 < x < -0.7, z = Zo, where the pulse is cho- 
sen so that it crosses a separatrix. For the single wave pulse, the 
boundary conditions were set such that: 

f . . riO;r / 1x1 r f -1 < < -0.7 

I otherwise 
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Fig. 8. Contours of numerical simulation of V for a fast wave sent in from lower boundary for -1 < x < l,z = 0.1 and its resultant 
propagation at times (a) ^=0.055, (b) t=0.n, (c) t=0A65, (d) t =0.22, (e) t=0.215 and (/) t=033, (g) ^=0.385, (h) t=0A4, (/) 
^=0.495, 0') t =0.55, (k) ^=0.605, (/) ^=0.66, (m) ^=0.72, (n) ^=0.74, (o) ^=0.76, (p) t =0.78, (q) r=0.8 and (r) ^=0.82, (s) t=OM, (t) 
t=OM, (u) r=0.88, (v) t =0.9, (w) ^=0.92 and (x) ^=0.94, labelling from top left to bottom right. 



d_ 

dx 



' — ^y 



d_ 

dx 



d 



0. 



z=1.6 



Also, zo = 0.2 was chosen (as opposed to zo =0.1) as it demon- 
strates the nature more clearly (zo = 0.1 gets a bit messy close to 
(x,z) = (-0.5,0)). 

We found that the linear Alfven wave propagates up from the 
z = Zo boundary and begins to spread out, following the field 
lines. This can be seen in Figure [121 As it follows the field lines, 
part of the wave accumulates along the separatrix and appears 
to stick to the separatrices (depositing parts of itself as it goes). 
This behaviour stretches the wave thinly along the separatrices. 



Part of the wave also appears to be propagating away from the 
magnetic skeleton (though it is actually just following the field 
lines and spreading out). These two goals for diff'erent parts of 
the wave result in a sharp wave front being formed. Note that 
as the wavefront gets very sharp (as seen in the bottom row of 
subfigures), the resolution of the simulation is unable to fully 
resolve the wave behaviour and accumulation. However, it is 
clear that part of the wave is propagating away from the mag- 
netic skeleton and the other part of the wave has accumulated 
along the separatrices (thus displaying the same nature as seen 
in Paper II). 
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Fig. 10. Analytical solution of V for WKB approximation of a fast wave sent in from lower boundary for -1 < x < 1, z = 0.1 and its 
resultant propagation at times (a) ^=0.055, (b) r=0.1 1, (c) t=0A65, (d) t =0.22, (e) t=0.215 and (/) t=033, (g) r=0.385, (h) t=0A4, 
(0 ^=0.495, (7) t =0.55, (k) r=0.605, (/) ^=0.66, (m) t=0J2, (n) ^=0.74, (o) r=0.76, (p) t =0.78, (q) r=0.8, (r) r=0.82, (s) ^=0.84, (0 
^=0.86, (u) ^=0.88, (v) t =0.9, (w) t=0.92 and (x) ^=0.94, labelling from top left to bottom right. The curves consisting of the crosses 
represent the front, middle and back edges of the WKB wave solution, where the pulse enters from the top of the box. 



4.2. WKB approximation 



Using the magnetic field for the two dipoles (equation|4l), we can 
apply the WKB approximation. Substituting V3; = ^''^^^'^^ • e~^^^ 
Again, we can use our WKB approximation to gain insight into into equation ([T3]) and assuming that » 1, leads to a first 
the numerical simulation. The Alfven equations we have to solve order PDE of the form (jc, z, 0, p, q) = 0. Applying the method 
are: of characteristics, we generate the equations: 



dfi 



d \ 



d 

' dx 



(13) ^ 
ds 
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Fig. 12. Contours of Vy for an Alfven wave sent in from lower boundary for -1 < x < -0.7 and its resultant propagation at times 
(a) r=0.05, (b) t=OA, (c) t=0A5, (d) t = 0.2, (e) t=OA, (/) ^=0.6, (g) t = 0.8, (h) t=0.9, (/) ^=1.0, (j) t = 1.1, (k) t=l2, (/) t=lA, 
labelling from top left to bottom right. 



ds \ dx dx I ' ds ^ ' 

dq _ IdB^ ^ ^^B^ ^ - b ^ (14) 

ds \ dz dz I ' ds ^ 

where co is the frequency of our wave, s is some parameter along 
the characteristic and ^ = {B^p + B^q). These five ODEs were 
solved numerically using a fourth-order Runge-Kutta method. 
By plotting positions of constant we can visualise the wave- 
front. This can be seen in Figure [131 Here, we have plotted the 
wave fronts at the same times as in the numerical simulation 
(Figure [121) to allow a direct comparison. The agreement is very 
good. Again, each wavefront is indicated by tiny crosses. Not 
only can we see part of the wave accumualting along the sepa- 
ratrices more clearly now, but we can also see how stretched the 
wave gets across the skeleton. 

We can also use our WKB solution to plot the particle paths of 
individual elements from the initial wave. These can be seen in 
Figure [T4l 

5. Conclusions 

This paper describes an investigation into the nature of MHD 
waves in the neighbourhood of two magnetic dipoles. From the 
work above, it has been seen that when a fast wave is generated 
on the lower boundary it propagates upwards and is influenced 
by the magnetic configuration. Firstly, it is influenced by the two 
regions of high Alfven speed; this distorts the initally straight 
wave to form two peaks (with maxima located over the loci of 
the magnetic field). These localised areas of high speed is a new 
eff'ect not seen previously in Papers I or II. 
The wave continues to propagate upwards but travels slower as 
it approaches the null (located at (x, z) = (0, 0.5)) due to the de- 
creasing Alfven speed. Also, the linear fast wave cannot cross 
the null (as its speed drops to zero at that point). Meanwhile, 
the rest of the wave (the wings) continue to propagate upwards 
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Fig. 14. Plots of WKB solution for an Alfven wave sent in from 
the lower boundary and its resultant particle paths (thick lines). 
Starting points of -1 < < -0.7 at intervals of ^ are plotted. 
Also shown in red is starting point = -0.758 



and spread out (the fast wave propagates roughly isotropically). 
Thus, this gives the visual impression that the wave is pivoting 
about the null point (which is acting as a fulcrum). As this pivot- 
ing continues, the two wings eventually cross each other and, due 
to the linear nature of the model, both wings pass each other (do 
not afl'ect each other). A full non-linear treatment of the equa- 
tions may reveal some very diff'erent behaviour. 
The two wings continue to pivot about their fulcrum. Close to 
the null, part of the wave is wrapping tightly around it. Near 
the null, a refraction efl'ect spirals the wave into the X-point, 
wrapping it around again and again. The rest of the wave (above 
the magnetic skeleton) appears to simply continue to rise and 
spread out. The wave is stretched between its two goals (part 
wrapping around the null and part rising away from the magnetic 
skeleton) and this leads to the wave splitting; this occurs near the 
regions of high Alfven speed (the localised high speed thins the 
wave and forces the split). Thus, part of the (now split) wave 
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Fig. 13. Analytical solution for WKB approximation of an Alfven wave sent in from lower boundary for -1 < x < -0.7, and its 
resultant propagation at times (a) t=0.05, (b) t=OA, (c) t=0A5, (d) t = 0.2, (e) t=OA, (/) ^=0.6, (g) t = 0.8, (h) t=0.9, (/) t=l.O, (j) 
t = 1.1, (^) ^= 1.2, (/) t=lA, labelling from top left to bottom right. The lines represent the front, middle and back edges of the WKB 
wave solution. 



spirals into the null and the other part propagates away from the 
magnetic skeleton; no longer heavily influenced by the magnetic 
null point. These parts form a distinctive cross shape. It was also 
demonstrated that 40% of the wave is trapped by the null, where 
the fraction captured is of the form j^x^^i^ (a, zo)- 

Thus, it has been shown that the key result for the fast wave 
from Paper I (i.e. that the fast wave accumulates at null points) 
persists in this more complication magnetic configuration, but 
only if the wave is close enough to the null point. This makes 
intuitive sense; if the fast wave is too far away from the mag- 
netic null, it will not feel its elfect (until it is much closer). The 
physical significance of this is that any fast magnetoacoustic dis- 
turbance near such a magnetic configuration will propagate and 
eventally split; the parts close enough to the null will be drawn 
towards it and focus all its energy at that point, the rest will prop- 
agate away from the system (taking its wave energy with it). The 
part trapped by the refraction eff'ect of the null will lead to large 
gradients in the perturbed magnetic field, and thus this is where 
current density build-up will occur. We believe ohmic dissipation 
will extract the energy in the wave at this point and wave heating 
will naturally occur at the null point (as in Paper I). However, 
this is only a claim as the numerical experiments and analytical 
work described above were all conducted in an ideal plasma. 

We also looked at the behaviour of the Alfven wave near our 
two dipoles configuration. We saw that when a linear Alfven 
wave was generated on the lower boundary, its propagation is 
confined to the path of the field line it starts on (the Alfven wave 
cannot cross field lines). This agrees with one of the key results 
found for Alfven waves in Paper II; the Alfven wave follows 
the field lines and each fluid element is trapped on the field line 
it starts on. Thus, part of the wave accumulates along the sep- 
aratrix (depositing parts of itself as it goes) and this behaviour 
stretches the wave thinly along the separatrices. The other part 
of the wave appears to be propagating away from the magnetic 
skeleton, though it is actually just following the field lines (and 



these are spreading out). This was apparent in the numerical sim- 
ulation but more clearly seen in the WKB approximation. These 
two goals for diff'erent parts of the wave result in a sharp wave- 
front being formed. Dissipation may extract the energy from the 
wave at this point. 

Furthermore, since we have the Alfven wave accumulating and 
stretched along the separatrices, we will have a build-up of cur- 
rent density at these locations. Therefore, since current is accu- 
mulating along the separatrices, this will be the location where 
ohmic dissipation will extract the energy in the wave. Thus, 
Alfven wave heating will naturally occur along the separatrices 
due to dissipation. The rest of the Alfven wave (the part trav- 
elling away from the magnetic skeleton) does not create large 
gradients and thus does not generate large currents. This part of 
the wave simply travels out of our system (taking its energy with 
it). 
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